Identification of small juvenile scombrids from northwest tropical Australia using mitochondrial DNA cytochrome b sequences

Small juveniles of the nine species of scombrids in Australian waters are morphologically similar to one another and, consequently, difficult to identify to species level. We show that the sequence of the mitochondrial DNA cytochrome b gene region is a powerful tool for identification of these young fish. Using this method, we identified 50 juvenile scombrids collected from Exmouth Bay, Western Australia. Six species of scombrids were apparent in this sample of fish: narrow-barred Spanish mackerel (Scomberomorus commerson), Indian mackerel (Rastrelliger kanagurta), frigate tuna (Auxis thazard), bullet tuna (Auxis rochei), leaping bonito (Cybiosarda elegans), and kawakawa (Euthynnus affinis). The presence of Indian mackerel, frigate tuna, leaping bonito, and kawakawa is the first indication that coastal waters may be an important spawning habitat for these species, although offshore spawning may also occur. The occurrence of small juvenile S. commerson was predicted from the known spawning patterns of that species, but other mackerel species (Scomberomorus munroi, Scomberomorus queenslandicus, Scomberomorus semifasiciatus) likely to be spawning during the sampling period were not detected among the 50 small juveniles analyzed here.

Mitochondrial DNA supports the identification of two endangered river sharks (Glyphis glyphis and Glyphis garricki) across northern Australia

The river sharks (genus Glyphis) are a small group of poorly known sharks occurring in tropical rivers and estuarine waters across northern Australia, south-east Asia and the subcontinent. The taxonomy of the genus has long been unclear due to very few individuals having been caught and examined, resulting in a paucity of data regarding their distribution, biology and ecology. Only recently has attention focussed on the two Australian species, G. glyphis and G. garricki. This study is a result of a rare opportunity to collate the few samples that have been collected from these species and the bull shark Carcharhinus leucas, which shares an overlapping range. These samples were analysed using the DNA barcoding approach (cox1 mitochondrial gene), compared with six other species of carcharhinids and evaluated in light of the current taxonomic classification. Nine species-specific nucleotide differences were found between G. glyphis and G. garricki and no intra-specific variation provides strong support for the separation into distinct species. Significant differences were also observed at the inter-generic level, with Glyphis forming a distinct clade from Carcharhinus. This study provides the basis for future molecular studies required to better address conservation issues confronting G. glyphis and G. garricki in Australia.

Mitochondrial DNA diversity of Cleruchoides noackae (Hymenoptera: Mymaridae): a potential biological control agent for Thaumastocoris peregrinus (Hemiptera: Thaumastocoridae)

Carpintero and Dellap, (Hemiptera: Thaumastocoridae) is a native Australian sap-feeding insect that has become invasive and seriously damaging to commercially grown in the Southern Hemisphere. Lin and Huber (Hymenoptera: Mymaridae) was recently discovered as an egg parasitoid of the Thaumastocoridae in Australia. Mitochondrial DNA (mtDNA; cytochrome oxidase subunit I, COI) sequence diversity amongst 104 individuals from these native populations revealed 24 sequence haplotypes. The COI haplotypes of individuals collected from the Sydney and Southeast Queensland clustered in distinct groups, indicating limited spread of the insect between the regions. Individuals collected from Perth in Western Australia were represented by four COI haplotypes. Although this population is geographically more isolated from other populations, two COI haplotypes were identical to haplotypes found in the Sydney region. The results suggest that has recently been introduced into Perth, possibly from the Sydney area. The high mtDNA diversity and limited spread that is suggested for is in contrast to the lack of geographic associated mtDNA diversity and extensive spread of . If implemented as a biological control agent, this factor will need to be considered in collecting and releasing .

Hybridisation, paternal leakage and mitochondrial DNA linearization in three anomalous fish (Scombridae)

Using mitochondrial DNA for species identification and population studies assumes that the genome is maternally inherited, circular, located in the cytoplasm and lacks recombination. This study explores the mitochondrial genomes of three anomalous mackerel. Complete mitochondrial genome sequencing plus nuclear microsatellite genotyping of these fish identified them as Scomberomorus munroi (spotted mackerel). Unlike normal S. munroi, these three fish also contained different linear, mitochondrial genomes of Scomberomorus semifasciatus (grey mackerel). The results are best explained by hybridisation, paternal leakage and mitochondrial DNA linearization. This unusual observation may provide an explanation for mtDNA outliers in animal population studies. © 2013.

DNA bar-coding reveals source and patterns of Thaumastocoris peregrinus invasions in South Africa and South America.

Thaumastocoris peregrinus is a recently introduced invertebrate pest of non-native Eucalyptus plantations in the Southern Hemisphere. It was first reported from South Africa in 2003 and in Argentina in 2005. Since then, populations have grown explosively and it has attained an almost ubiquitous distribution over several regions in South Africa on 26 Eucalyptus species. Here we address three key questions regarding this invasion, namely whether only one species has been introduced, whether there were single or multiple introductions into South Africa and South America and what the source of the introduction might have been. To answer these questions, bar-coding using mitochondrial DNA (COI) sequence diversity was used to characterise the populations of this insect from Australia, Argentina, Brazil, South Africa and Uruguay. Analyses revealed three cryptic species in Australia, of which only T. peregrinus is represented in South Africa and South America. Thaumastocoris peregrinus populations contained eight haplotypes, with a pairwise nucleotide distance of 0.2-0.9% from seventeen locations in Australia. Three of these haplotypes are shared with populations in South America and South Africa, but the latter regions do not share haplotypes. These data, together with the current distribution of the haplotypes and the known direction of original spread in these regions, suggest that at least three distinct introductions of the insect occurred in South Africa and South America before 2005. The two most common haplotypes in Sydney, one of which was also found in Brisbane, are shared with the non-native regions. Sydney populations of T. peregrinus, which have regularly reached outbreak levels in recent years, might thus have served as source of these three distinct introductions into other regions of the Southern Hemisphere.

The complete validated mitochondrial genome of the silver gemfish Rexea solandri (Cuvier, 1832) (Perciformes, Gempylidae)

The silver gemfish Rexea solandri is an important economic resource but vulnerable to overfishing in Australian waters. The complete mitochondrial genome sequence is described from 1.6 million reads obtained via next generation sequencing. The total length of the mitogenome is 16,350 bp comprising 2 rRNA, 13 protein-coding genes, 22 tRNA and 2 non-coding regions. The mitogenome sequence was validated against sequences of PCR fragments and BLAST queries of Genbank. Gene order was equivalent to that found in marine fishes.

The complete validated mitochondrial genome of the black marlin Istiompax indica (Cuvier, 1832)

Two complete mitochondrial genomes of the black marlin Istiompax indica were assembled from approximately 3.5 and 2.5 million reads produced by Ion Torrent next generation sequencing. The complete genomes were 16,531 bp and 16,532 bp in length consisting of 2 rRNA, 13 protein-coding genes, 22tRNA and 2 coding regions. They demonstrated a similar A + T base (52.6%) to other teleosts. Intraspecific sequence variation was 99.5% for three I. indica mitogenomes and 99.7% for X. gladius. A lower value (85%) was found for the I. platypterus mitogenomes from genbank and accredited to inadvertent inclusion of gene regions from a con-familial species in one record, highlighting the need for cautious downstream use of genbank data. © 2014 Informa UK Ltd.

A mitochondrial species identification assay for Australian blacktip sharks (Carcharhinus tilstoni, C. limbatus and C. amblyrhynchoides) using real-time PCR and high-resolution melt analysis

Tropical Australian shark fisheries target two morphologically indistinguishable blacktip sharks, the Australian blacktip (Carcharhinus tilstoni) and the common blacktip (C. limbatus). Their relative contributions to northern and eastern Australian coastal fisheries are unclear because of species identification difficulties. The two species differ in their number of precaudal vertebrae, which is difficult and time consuming to obtain in the field. But, the two species can be distinguished genetically with diagnostic mutations in their mitochondrial DNA ND4 gene. A third closely related sister species, the graceful shark C. amblyrhynchoides, can also be distinguished by species-specific mutations in this gene. DNA sequencing is an effective diagnostic tool, but is relatively expensive and time consuming. In contrast, real-time high-resolution melt (HRM) PCR assays are rapid and relatively inexpensive. These assays amplify regions of DNA with species-specific genetic mutations that result in PCR products with unique melt profiles. A real-time HRM PCR species-diagnostic assay (RT-HRM-PCR) has been developed based on the mtDNA ND4 gene for rapid typing of C. tilstoni, C. limbatus and C. amblyrhynchoides. The assay was developed using ND4 sequences from 66 C. tilstoni, 33. C. limbatus and five C. amblyrhynchoides collected from Indonesia and Australian states and territories; Western Australia, the Northern Territory, Queensland and New South Wales. The assay was shown to be 100% accurate on 160 unknown blacktip shark tissue samples by full mtDNA ND4 sequencing.

Phylogeography of the green turtle, Chelonia mydas, in the Southwest Indian Ocean

Patterns of mitochondrial DNA (mtDNA) variation were used to analyse the population genetic structure of southwestern Indian Ocean green turtle (Chelonia mydas) populations. Analysis of sequence variation over 396 bp of the mtDNA control region revealed seven haplotypes among 288 individuals from 10 nesting sites in the Southwest Indian Ocean. This is the first time that Atlantic Ocean haplotypes have been recorded among any Indo-Pacific nesting populations. Previous studies indicated that the Cape of Good Hope was a major biogeographical barrier between the Atlantic and Indian Oceans because evidence for gene flow in the last 1.5 million years has yet to emerge. This study, by sampling localities adjacent to this barrier, demonstrates that recent gene flow has occurred from the Atlantic Ocean into the Indian Ocean via the Cape of Good Hope. We also found compelling genetic evidence that green turtles nesting at the rookeries of the South Mozambique Channel (SMC) and those nesting in the North Mozambique Channel (NMC) belong to separate genetic stocks. Furthermore, the SMC could be subdivided in two different genetic stocks, one in Europa and the other one in Juan de Nova. We suggest that this particular genetic pattern along the Mozambique Channel is attributable to a recent colonization from the Atlantic Ocean and is maintained by oceanic conditions in the northern and southern Mozambique Channel that influence early stages in the green turtle life cycle.

A universal test to determine the integrity of RNA, and its suitability for reverse transcription, in animal tissue laboratory specimens

Degradation of RNA in diagnostic specimens can cause false-negative test results and potential misdiagnosis when tests rely on the detection of specific RNA sequence. Current molecular methods of checking RNA integrity tend to be host species or group specific, necessitating libraries of primers and reaction conditions. The objective here was to develop a universal (multi-species) quality assurance tool for determining the integrity of RNA in animal tissues submitted to a laboratory for analyses. Ribosomal RNA (16S rRNA) transcribed from the mitochondrial 16S rDNA was used as template material for reverse transcription to cDNA and was amplified using polymerase chain reaction (PCR). As mitochondrial DNA has a high level of conservation, the primers used were shown to reverse transcribe and amplify RNA from every animal species tested. Deliberate degradation of rRNA template through temperature abuse of samples resulted in no reverse transcription and amplification. Samples spiked with viruses showed that single-stranded viral RNA and rRNA in the same sample degraded at similar rates, hence reverse transcription and PCR amplification of 16S rRNA could be used as a test of sample integrity and suitability for analysis that required the sample's RNA, including viral RNA. This test will be an invaluable quality assurance tool for determination of RNA integrity from tissue samples, thus avoiding erroneous test results that might occur if degraded target RNA is used unknowingly as template material for reverse transcription and subsequent PCR amplification.

The genetic structure of Australasian green turtles (Chelonia mydas): Exploring the geographical scale of genetic exchange

Ecological and genetic studies of marine turtles generally support the hypothesis of natal homing, but leave open the question of the geographical scale of genetic exchange and the capacity of turtles to shift breeding sites. Here we combine analyses of mitochondrial DNA (mtDNA) variation and recapture data to assess the geographical scale of individual breeding populations and the distribution of such populations through Australasia. We conducted multiscale assessments of mtDNA variation among 714 samples from 27 green turtle rookeries and of adult female dispersal among nesting sites in eastern Australia. Many of these rookeries are on shelves that were flooded by rising sea levels less than 10 000 years (c. 450 generations) ago. Analyses of sequence variation among the mtDNA control region revealed 25 haplotypes, and their frequency distributions indicated 17 genetically distinct breeding stocks (Management Units) consisting either of individual rookeries or groups of rookeries in general that are separated by more than 500 km. The population structure inferred from mtDNA was consistent with the scale of movements observed in long-term mark-recapture studies of east Australian rookeries. Phylogenetic analysis of the haplotypes revealed five clades with significant partitioning of sequence diversity (Φ = 68.4) between Pacific Ocean and Southeast Asian/Indian Ocean rookeries. Isolation by distance was indicated for rookeries separated by up to 2000 km but explained only 12% of the genetic structure. The emerging general picture is one of dynamic population structure influenced by the capacity of females to relocate among proximal breeding sites, although this may be conditional on large population sizes as existed historically across this region.

Genetic population structure of red snappers (Lutjanus malabaricus Bloch & Schneider, 1801 and Lutjanus erythropterus Bloch, 1790) in central and eastern Indonesia and northern Australia

The genetic population structure of red snapper Lutjanus malabaricus and Lutjanus erythropterus in eastern Indonesia and northern Australia was investigated by allozyme electrophoresis and sequence variation in the control region of mtDNA. Samples were collected from eight sites in Indonesia and four sites in northern Australia for both species. A total of 13 allozyme loci were scored. More variable loci were observed in L. malabaricus than in L. erythropterus. Sequence variation in the control region (left domain) of the mitochondrial genome was assessed by RFLP and direct sequencing. MtDNA haplotype diversity was high (L. erythropterus, 0.95 and L. malabaricus, 0.97), as was intraspecific sequence divergence, (L. erythropterus, 0.0-12.5% and L. malabaricus, 0.0-9.5%). The pattern of mtDNA haplotype frequencies grouped both species into two broad fisheries stocks with a genetic boundary either between Kupang and Sape (L. malabaricus) or between Kupang and Australian Timor Sea (L. erythropertus). The allozyme analyses revealed similar boundaries for L. erythropterus. Seven allozymes stocks compared to two mtDNA stocks of L. malabaricus including Ambon, which was not sampled with mtDNA, however, were reported. Possible reasons for differences in discrimination between the methods include: i) increased power of multiple allozyme loci over the single mtDNA locus, ii) insufficient gene sampling in the mtDNA control region and iii) relative evolutionary dynamics of nuclear (allozyme loci) and mitochondrial DNA in these taxa. Allozyme and haplotype data did not distinguish separate stocks among the four Australian locations nor the central Indonesian (Bali and Sape locations) for both L. malabaricus and L. erythropterus.

Concordant molecular and phenotypic data delineate new taxonomy and conservation priorities for the endangered mountain yellow-legged frog

The mountain yellow-legged frog Rana muscosa sensu lato, once abundant in the Sierra Nevada of California and Nevada, and the disjunct Transverse Ranges of southern California, has declined precipitously throughout its range, even though most of its habitat is protected. The species is now extinct in Nevada and reduced to tiny remnants in southern California, where as a distinct population segment, it is classified as Endangered. Introduced predators (trout), air pollution and an infectious disease (chytridiomycosis) threaten remaining populations. A Bayesian analysis of 1901 base pairs of mitochondrial DNA confirms the presence of two deeply divergent clades that come into near contact in the Sierra Nevada. Morphological studies of museum specimens and analysis of acoustic data show that the two major mtDNA clades are readily differentiated phenotypically. Accordingly, we recognize two species, Rana sierrae, in the northern and central Sierra Nevada, and R. muscosa, in the southern Sierra Nevada and southern California. Existing data indicate no range overlap. These results have important implications for the conservation of these two species as they illuminate a profound mismatch between the current delineation of the distinct population segments (southern California vs. Sierra Nevada) and actual species boundaries. For example, our study finds that remnant populations of R. muscosa exist in both the southern Sierra Nevada and the mountains of southern California, which may broaden options for management. In addition, despite the fact that only the southern California populations are listed as Endangered, surveys conducted since 1995 at 225 historic (1899-1994) localities from museum collections show that 93.3% (n=146) of R. sierrae populations and 95.2% (n=79) of R. muscosa populations are extinct. Evidence presented here underscores the need for revision of protected population status to include both species throughout their ranges.

Grey Mackerel (Scomberomorus semifasciatus) microsatellite genotypes (nine loci) representing twelve regions along Australia's northern coastline from the eastern coast of Queensland to the western coast of Western Australia.

Scomberomorus semifasciatus is an Australian endemic found in tropical, coastal waters from eastern to western Australia. Commercial and recreational exploitation is common and regulated by state-based authorities. This study used mitochondrial DNA sequence and microsatellite markers to elucidate the population structure of Scomberomorus semifasciatus collected from twelve, equidistant sampling locations. Samples (n=544) were genotyped with nine microsatellite loci, and 353 were sequenced for d-loop (384 bp) and ATP (800bp) mitochondrial DNA gene regions. Combined interpretation of microsatellite and mtDNA data identified four genetic stocks of S. semifasciatus: Western Australia, northwest coast of the Northern Territory, Gulf of Carpentaria and the east coast of Queensland. Connectivity among stocks across northern Australia from the Northern Territory to the east coast of Queensland was high, but in contrast, there was a clear genetic break between populations in Western Australia compared to the rest of northern Australia. This indicates a restriction to gene flow possibly associated with suboptimal habitat along the Kimberley coast (northwestern Australia). The appropriate scale of management for this species corresponds to the jurisdictions of the three Australian states, except that the Gulf of Carpentaria stock should be co-managed by authorities in Queensland and Northern Territory.

Migration of green turtles (Chelonia mydas) from Australasian feeding grounds inferred from genetic analyses

Coastal seagrass habitats in tropical and subtropical regions support aggregations of resident green turtles (Chelonia mydas) from several genetically distinct breeding populations. Migration of individuals to their respective dispersed breeding sites provides a complex pattern of migratory connectivity among nesting and feeding habitats of this species. An understanding of this pattern is important in regions where the persistence of populations is under threat from anthropogenic impacts. The present study uses mitochondrial DNA and mixed-stock analyses to assess the connectivity among seven feeding grounds across the north Australian coast and adjacent areas and 17 genetically distinct breeding populations from the Indo-Pacific region. It was hypothesised that large and geographically proximate breeding populations would dominate at nearby feeding grounds. As expected, each sampled feeding area appears to support multiple breeding populations, with two aggregations dominated by a local breeding population. Geographic distance between breeding and feeding habitat strongly influenced whether a breeding population contributed to a feeding ground (wi = 0.654); however, neither distance nor size of a breeding population was a good predictor of the extent of their contribution. The differential proportional contributions suggest the impact of anthropogenic mortality at feeding grounds should be assessed on a case-by-case basis.